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In 1999, the highly compact and variable BL Lac object AO 0235+164 was identified as the highest 
brightness temperature active galactic nucleus observed with the VLBI Space Observatory Programme 
(VSOP), with Tb > 5.8 x 10 13 K (Frey et al. 2000). The sub-milliarcsecond radio structure of this source 
ly*-) \ has been studied with dual- frequency (1.6 and 5 GHz), polarization-sensitive VSOP observations during 

2001 and 2002. Here we present the results of this monitoring campaign. At the time of these observations, 
the source was weakly polarized and characterized by a radio core that is clearly resolved on space-ground 
baselines. 

Key words: BL Lacertae objects: individual (AO 0235+164) — galaxies: active — radio continuum: 
galaxies — techniques: interferometric 
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■ 1. Introduction et al. (2001) recently claimed to identify two components 

. £h ! well within 1 mas of the core that show apparent super- 

^ . AO 0235+164 is an extensively studied BL Lac ob- luminal motion with speeds up to /3 a pp = 30/i 1 c. This 
ject at redshift z = 0.940 (Cohen et al. 1987). The ob- would imply a Doppler factor of at least 90 in the jet when 
ject is seen through a foreground group of galaxies at it points directly to the observer. Such high Doppler fac- 
z = 0.524 responsible for spectral lines in both absorption tors (~ 100) were reported by Kraus et al. (1999) based 
and emission. Absorption lines from z = 0.851 are also on radio variability measurements, and by Fujisawa et 
present in the spectrum. The source is highly variable al. (1999) derived from VLBI and radio total flux density 
over the whole electromagnetic spectrum from radio to observations. The result is also consistent with our earlier 
7-rays. Variability time scales from years to days are ob- first-epoch 5-GHz VLBI Space Observatory Programme 
served, as well as intra-day variability. See e.g. Raiteri et (VSOP) imaging observation (Frey et al. 2000). We de- 
al. (2001) and references therein for a detailed account of rived a lower limit to the brightness temperature, Tb > 
the optical and radio variability studies of AO 0235+164 5.8 x 10 13 K, based on an unresolved core component, 
over the past three decades. This is the highest value measured directly and implies 

The high-resolution radio structure of the source at a Doppler factor of ~ 100. 
milli-arcsecond (mas) angular scales revealed by ground- However, the Doppler boosting also varies with time as 
based Very Long Baseline Interferometry (VLBI) is char- can be inferred from the brightness temperatures derived 
acterized by a dominant compact core. Occasionally, faint from VLBI data (Frey et al. 2000). This, and the large 
extensions are seen at various position angles mainly be- changes in the apparent VLBI jet position angle are qual- 
tween the north and east. A collection of references to itativcly well explained by small deviations in a jet that 
earlier VLBI observations is given by Frey et al. (2000). intrinsically points very close to the line of sight. Indeed, 
Based on a series of 43-GHz Very Long Baseline Array Ostorero et al. (2004) interpreted the long-term quasi- 
(VLBA) observations spanning almost two years, Jorstad periodic variability of AO 0235+164 at multiple wave- 
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Table 1. Space VLBI observations of AO 0235+164 at 
1.6 GHz (top) and 5 GHz (bottom). Polarization-sensitive 
experiments are marked with asterisks 





Epoch 


Ground network 


Correlator 


A 


1999 Jan 31 


AT HH NO SH (4) 


Penticton 


B* 


2001 Feb 4 


VLBA AR (11) 


Socorro 


c* 


2001 Aug 3 


VLBA AR GO (12) 


Socorro 


D 


1999 Feb 1 


AT SH UD (3) 


Penticton 


E* 


2001 Feb 2 


VLBA Y (11) 


Socorro 


p* 


2001 Aug 2 


VLBA AR (11) 


Socorro 


G 


2001 Aug 14 


AR AT UD (3) 


Mitaka 


H* 


2002 Jan 26 


VLBA (9) 


Socorro 



Antenna names are AR: Arecibo (Puerto Rico); AT: 
Australia Telescope Compact Array; GO: Goldstone 
(USA); HH: Hartebeesthoek (South Africa); NO: Noto 
(Italy); SH: Shcshan (China); UD: Usuda (Japan); 
VLBA: Very Long Baseline Array (USA); Y: Very Large 
Array (VLA, USA) 

bands as orientation variations in a helical jet. In their 
model, non-periodic features are explained by flow in- 
stabilities. Short timcscale (intra-hour) variability data 
taken in 2000 arc best interpreted by invoking interstellar 
scintillation of a source of at least 0.015 mas in angular 
size, implying a Doppler factor smaller than 40 in the rel- 
ativistically moving component (Peng & dc Bruyn 2004). 

As opposed to other BL Lac objects, high-resolution 
radio polarization studies using ground-based VLBI at 
5 GHz (Gabuzda et al. 1992) and 8.4 GHz (Gabuzda & 
Cawthorne 1996) found no convincing structure in either 
total or polarized intensity apart from the compact core of 
AO 0235+164. Our VSOP monitoring observations pre- 
sented here were aimed at investigating how the sub-mas 
scale radio structure of AO 0235+164 varies, and to see 
whether linear polarization structure can be detected with 
the superior angular resolution offered by VSOP at the fre- 
quencies of 1.6 and 5 GHz. The VSOP satellite, HALCA, 
receives only left-circularly polarized radiation, and has a 
limited sensitivity compared to ground-based VLBI sta- 
tions. However, dual-polarization experiments with the 
co-observing ground network (e.g. the VLBA) allow po- 
larization imaging of sources that show sufficiently high 
correlated polarized flux density (Kemball et al. 2000). 

2. Observations and data reduction 

The details of the 3-year Space VLBI (SVLBI) monitor- 
ing program are summarized in Table 1. AO 0235+164 
was observed at 1.6 and 5 GHz with VSOP at a to- 
tal of 8 epochs, including the two experiments in 1999 
(marked with A and D) already published by Frey et 
al. (2000) . Here we focus on the results of the other obser- 
vations. For each experiment, the data from the space ra- 
dio telescope on board the HALCA satellite (Hirabayashi 
et al. 2000a) were recorded at a subset of the tracking sta- 
tions at Goldstone, Green Bank (USA), Robledo (Spain), 
Tidbinbilla (Australia) and Usuda (Japan). The band- 


















CM 


- 




°o 

% o 
o 

° 

° % 


o 
o 




o 

- 

ft 








i 

% o o 


oo 
? °o 

8 


o c$S 


^ n od& o 

PL » „° o? 

r % o ° oPiP- 

o 








<*> 




»»o? 


O 

O 








O CP 


o 


8°° 










6° % 
« 




s 
























t 




t t 


t 




I 




D 




E F 


H 




Q 


O 

























1996 1998 2000 2002 2004 

time (year) 

Fig. 1. Total flux density history of AO 0235+164 at 4.8 GHz 
between 1996 and 2005. Two- week average UMRAO measure- 
ments are shown. Our SVLBI monitoring epochs are marked. 

width was 32 MHz. As shown in Table 1, on different occa- 
sions, different VSOP data processors in Mitaka (Japan), 
Penticton (Canada) and Socorro (USA) were used. The 
total number of ground stations for each experiment is 
given in parentheses in Table 1. Polarization-sensitive ob- 
servations involving 10 (in one case 9) antennas of the 
VLBA are marked with asterisks after the experiment 
identifier in the first column. 

Note that AO 0235+164 is included in the VSOP 
5-GHz Active Galactic Nucleus (AGN) Survey sample 
(Hirabayashi et al. 2000b). However, due to its variability, 
at the time of the 5-GHz prc-launch VLBA observations 
(Fomalont et al. 2000) it failed to reach the required cor- 
related flux density for inclusion in the final source list to 
be observed with HALCA in the Survey. 

At each observing epoch, multi-frequency measure- 
ments of the total radio flux density were obtained quasi- 
contemporaneously at the University of Michigan Radio 
Astronomy Observatory (UMRAO) and the Arecibo 
Observatory. AO 0235+164 is included in the long-term 
UMRAO monitoring program of total and polarized flux 
densities at 4.8, 8 and 14.5 GHz (Aller et al. 1999; Aller et 
al. 2005). An up-to-date 4.8-GHz total flux density curve 
covering our monitoring period is given in Fig. 1. See 
Aller et al. (2005) for multiple frequencies and a longer 
time range. The single-dish flux density measurements 
were invaluable for verifying the calibration of our inter- 
ferometer data. 

The US National Radio Astronomy Observatory 
(NRAO) Astronomical Image Processing System (AIPS) 
was used for the data calibration and imaging. The visibil- 
ity amplitudes were calibrated using system temperatures 
measured at the antennas wherever available. Nominal 
values were applied for HALCA and occasionally for some 
ground antennas (AR, AT and UD). Due to the sufficiently 
large ground network in most experiments, these could be 
verified and adjusted in the calibration procedure. The 
initial phase calibration was done "manually" using short 
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Fig. 2. 1.6-GHz SVLBI images of AO 0235+164 taken on 
2001 Feb 4 and 2001 Aug 3 (epochs B and C). The first total 
intensity contours arc drawn at ±2 mjy/bcam in each image. 
The positive contour levels increase by a factor of 2. The peak 
brightnesses are 703 mjy/bcam (B) and 753 mjy/bcam (C). 
The restoring beams are 6.40 mas X 3.07 mas at the position 
angle (PA) -1° (B) and 5.20 mas X 2.00 mas at PA=5° (C), 
as displayed in the lower-left corners. The off-source RMS 
image noise (ltr) is 0.4 mjy/beam (B) and 0.5 mjy/beam 
(C). The right ascension (horizontal axis) and declination 
(vertical axis) scales have the brightness peak as origin. 



scans of data. 

In the case of the polarization-sensitive experiments (B, 
C, E, F and H), the right-left delay difference was re- 
moved with the AIPS procedure VLBACPOL. Fringe- 
fitting was performed with the AIPS task FRING over 5- 
minutc solution intervals. The total intensity images are 
the results of a hybrid mapping procedure involving sev- 
eral iterations with the tasks CALIB and IMAGR. The 
fourth root of the visibility data weights were used in a 
uniform weighting scheme in order to increase the relative 
importance of the baselines to the orbiting radio telescope 
and hence to improve the angular resolution. The total 
intensity images arc displayed in Fig. 2 (1.6 GHz) and 
Fig. 3 (5 GHz). Identical coordinate scales and contour 
levels are used for each epoch in both figures. No image 
is presented for the 5-GHz observation on 2001 Aug 14 
(experiment G), because only single baselines were avail- 
able for most of the time. However, the source was clearly 
detected on each baseline and the data are consistent with 
a slightly resolved compact component. 

For the polarization experiments, the task LPCAL was 
used to solve for instrumental polarization terms for the 
antennas. Absolute calibration of the polarization posi- 
tion angle was possible for experiments E and F. The cal- 
ibrator source 0048-097 was observed with the ground an- 
tennas. Integrated polarization measurements for the cal- 
ibrator were obtained from the NRAO Very Large Array 
(VLA) involved in the SVLBI experiment as a phased ar- 
ray (epoch E), and from the UMRAO monitoring (epoch 
F). In the latter case, due to a 2-month period between 
our epoch and the data base entry, the calibration is more 
uncertain. The corresponding electric vectors are shown 
in Fig. 3. The electric vector position angles are —11° and 
34° for experiments E and F, respectively. 

While most of the ground antennas (the VLBA. the 
phased VLA) recorded in both right (R) and left (L) 
circular polarization, HALCA and some of the ground- 
based antennas received only left circular polarization. 
Therefore only LL and RL correlations were obtained for 
certain baselines including those to the orbiting antenna. 
In AIPS, the procedure CXPOLN and the task CXCLN 
were used to perform complex polarization imaging. 



Results and discussion 



At 1.6 GHz, the total intensity images of AO 0235+164 
(Fig. 2) obtained at epochs B and C separated by 6 months 
in time show a resolved and somewhat extended structure. 
The low-brightness extension is on the eastern side of the 
compact core. The fractional polarization of the core was 
<1%- 

At 5 GHz, in contrast to our earlier results (epoch 
D, Frey et al. 2000), the source is well resolved on the 
space-ground baselines at three different epochs. This is 
clearly indicated by the plots of the correlated flux den- 
sity as a function of the projected baseline length (Fig. 4). 
Although there are hints of sub-mas scale structure in the 
images, these arc difficult to interpret due to the highly 
elongated restoring beam, a consequence of the asymmet- 



4 



S. Frey et al. 





H 




£ 
S 

3 0.8 



S 

% 0.8 



1.0 

| 

| 0.8 



100 200 300 400 

Mega Wavlngth 



500 600 



100 200 300 400 500 600 

Mega Wavlngth 



H 



Milli ARC SEC 



300 400 
Mega Wavlngth 



Fig. 3. 5-GHz SVLBI images of AO 0235+164 taken on 2001 
Feb 2, 2001 Aug 2 and 2002 Jan 26 (epochs E, F and H). The 
first total intensity contours arc drawn at ±4 mjy/beam in 
each image. The positive contour levels increase by a fac- 
tor of 2. The peak brightnesses are 859 mjy/beam (E), 
1243 mjy/beam (F) and 1008 mjy/beam (H). The restoring 
beams are 1.61 mas X 0.25 mas at PA=— 34° (E), 1.82 mas 
X 0.63 mas at PA=8° (F) and 2.23 mas X 0.27 mas at 
PA=53° (H), as displayed in the lower-left corners. The 
off-source RMS image noise (lc) is 0.9 mjy/beam (E, F) and 
1.1 mjy/beam (H). Electric vectors are superimposed (E and 
F), 1 mas corresponds to 10 mjy/beam polarized intensity. 



Fig. 4. Self-calibrated 5-GHz correlated flux density vs. pro- 
jected baseline length for the SVLBI observations of AO 
0235+164 taken on 2001 Feb 2, 2001 Aug 2 and 2002 Jan 
26 (epochs E, F and H). Only every tenth point is plot- 
ted for clarity. Ground-only baselines typically extend up 
to ~ 150 million wavelengths (MA). Total (single-dish) 
flux density values obtained at the UMRAO at nearly 
the same time are indicated with horizontal dashed lines. 
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Fig. 5. 5-GHz image of AO 0235+164 made from the 
VLA data only (2001 Feb 2). The first total inten- 
sity contours are drawn at ±0.25 mjy/beam. The posi- 
tive contour levels increase by a factor of 2. The peak 
brightness is 1249 mjy/beam. The restoring beam is 
1.25" X 0.71" at PA=-89°. The off-source RMS image 
noise (lcr) is 0.05 mjy/beam. Electric vectors are superim- 
posed, 1" corresponds to 10 mjy/beam polarized intensity. 

ric (u,w)-coverage. In fact, the source brightness distribu- 
tion is well fitted with a single circular Gaussian compo- 
nent of 0.20, 0.30 and 0.23 mas (FWHM) in the case of 
experiments E, F and H, respectively. These values imply 
brightness temperatures of ~ 10 12 K, more than an order 
of magnitude less than that derived in 1999, shortly af- 
ter a major total flux density outburst (Frey et al. 2000). 
The source has apparently become extended and, accord- 
ing to the UMRAO total flux density monitoring data, 
also faded since 1999. Our monitoring epochs fall into a 
long, relatively quiescent period in the flux density history 
of AO 0235+164 (Fig. 1). 

The nearly contemporaneous 5-GHz total flux density 
measurements are indicated with horizontal dashed lines 
in the correlated flux density plots (Fig. 4) for compari- 
son. It is interesting to note that the correlated flux densi- 
ties on the shortest VLBI baselines are considerably lower 
than the total flux densities, by up to ~ 150 mJy (epochs 
E, F). This phenomenon is often encountered at different 
epochs and frequencies by other authors as well. In partic- 
ular, Ros (2004) lists AO 0235+164 as one of the "compact 
but slightly resolved" sources found in the VLBA 2-cm 
Survey sample, suggesting that a weak extended emission 
is associated with a jet structure seen more or less face- 
on. In the framework of this survey, Kovalev et al. (2005) 
derived T B = 9.08 x 10 11 K at 15 GHz at the epoch of 2001 
Mar 15, consistently with our results discussed here. 

Since the phased VLA was included in the ground net- 
work in our 5-GHz experiment on 2001 Feb 2 (epoch E), 
we were also able to image AO 0235+164 using only VLA 
data. Our VLA image (Fig. 5) with about 3 orders of 
magnitude less resolution than that of the SVLBI im- 
ages shows weak extended emission surrounding the core, 
reaching 0.1% of the peak brightness in a component at 



an angular distance of 7"5 to the north- northwest. The 
structure seen is similar to earlier VLA results (Murphy 
et al. 1993). Unlike the VLBI data, the VLA observations 
fully account for the total flux density. About 10 mJy flux 
density can be attributed to the extended structure in the 
VLA image. We believe that diffuse low surface brightness 
emission is present on intermediate angular scales, but is 
resolved out by VLBI. 

As argued by e.g. Fujisawa et al. (1999), the radio jet 
in AO 0235+164 may lie very close (within 1°) to the line 
of sight. It is possible that this angle changes with a pe- 
riod of 5 — 6 years (Raitcri et al. 2001), as in the helical 
jet model proposed by Ostorero et al. (2004). However, 
the next major radio outburst predicted by the model 
for around early 2004 has yet to be detected (Raiteri et 
al. 2005). Qualitatively, this explains a wide range of phe- 
nomena observed in the source, including (some) variabil- 
ity, rapid changes in the VLBI jet component position 
angles (when observed), and the diffuse radio emission 
surrounding the core at larger scales. In this picture, our 
SVLBI monitoring covers a period post-1999 that cor- 
responds to a larger jet angle with respect to the line 
of sight, and therefore a less dominant Doppler boost- 
ing. With the jet Doppler factor S ~ 100 determined in 
1999, the Lorentz factor characterising the bulk plasma 
motion is 7 ^ 50, and the corresponding jet viewing angle 
is 9 < 0?5 (Urry & Padovani 1995, Appendix A). At our 
later epochs, the Doppler boosting decreased by a factor of 
~ 20. Under the simple assumption that the bulk Lorentz 
factor remained similar, the observed change in the bright- 
ness temperature can be explained by an increase in as 
small as < 5° . 

At our monitoring epochs, the fractional polarization 
of AO 0235+164 was quite low, ~ 0.5 - 2%. This result 
is fully consistent with the corresponding UMRAO mon- 
itorig data, where the polarized flux density was mea- 
sured to be around the detection limit at these epochs. 
Historically, values as high as ~ 5% have been measured 
for AO 0235+164. Besides the low level of polariza- 
tion, the structure found is simple, no feature apart from 
the compact core is being seen in our polarization im- 
ages. This rather atypical property of AO 0235+164 com- 
pared to other BL Lac objects has already been noted by 
Gabuzda et al. (1992) and Gabuzda & Cawthorne (1996) 
using ground-based VLBI data. In this "inactive" phase 
the degree of ordering of the magnetic field in the source 
is low and no propagating shock is emerging in the plasma 
flow. 

SVLBI observations of other BL Lac objects (Pushkarev 
et al. 2005, and references therein) resolve the polariza- 
tion structure of the inner jet components. This is typ- 
ically aligned with the jet axis, suggesting that the jets 
are associated with helical magnetic fields that propogate 
outward with the jet flow. 
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